Characterizing Genuine Multisite Entanglement in Isotropic Spin Lattices 
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We consider a class of large superposed states, obtained from dimer coverings on spin- 1/2 isotropic 
lattices. We show that they are genuinely multiparty entangled, irrespective of the geometry and 
dimension of the isotropic lattice. We then present an efficient method to characterize the genuine 
multisite entanglement in the case of isotropic square spin-1/2 lattices, with short-range dimcr 
coverings. We use this iterative analytical method to calculate the multisite entanglement of finite- 
sized lattices, which, through finite-size scaling, enables us to obtain the estimate of the multisite 
entanglement of the infinite square lattice. The method can be a useful tool to investigate other 
single- and multisite properties of such states. 



I. INTRODUCTION 

Genuine multisite entanglement is an important re- 
source in quantum information protocols and is known to 
offer significant advantage in quantum tasks in compari- 
son to bipartite entanglement [I]. In particular, it is the 
basic ingredient in measurement-based quantum compu- 
tation [2], and is beneficial in various quantum commu- 
nication protocols [3], including secret sharing [4] (cf. 
[5]). Apart from the conventional information tasks, the 
study of multisite entanglement turns out to be impor- 
tant in understanding many-body phenomena like quan- 
tum phase transitions [6] and in biological mechanisms 
to understand transport properties in the evolution of 
photsynthctic complexes [7]. Although bipartite entan- 
glement in the case of two spin-1/2 particles is rather 
well-understood, the situation is quite different in the 
case of classification and quantification of entanglement 
in higher dimensions as well as in multiparty systems. 
The fact that many-particle systems can have different 
types of useful entanglement, depending on the particu- 
lar information processing protocols under study, makes 
the quantification a formidable task. While there are 
several multiparty entanglement measures known in the 
literature [1, 8], it is in general difficult to compute them. 
However, for pure multisite states, it is possible to use the 
generalized geometric measure [9] , which is a computable 
measure of genuine multisite entanglement. 

The isotropic spin-1/2 lattice states, with short-range 
dimer coverings, are of considerable interest in investi- 
gating high-temperature superconductivity [10, 11], and 
critical phenomena in many body systems [12-14]. These 
short-range dimcr states can be efficiently simulated in 
laboratories using atoms in optical lattices [15] or in cav- 
ities using interacting photons [10]. The large number 
of superpositions of the nearest neighbor dimer coverings 
even in moderately large two-dimensional (2D) lattices, 
makes the computation of both bipartite and multisite 
entanglement of such spin states an arduous task. 

In this article, we consider quantum states, with dimer 
coverings, of an isotropic spin-1/2 lattice. We prove 
that such spin-1/2 rotationally invariant states are al- 
ways genuinely multisite entangled, irrespective of the 



lattice geometry and dimension. Wc present a method 
to analytically calculate the genuine multipartite entan- 
glement, viz. the generalized geometric measure [9], of 
these quantum spin states, with short-ranged dimer cov- 
erings, on a square lattice with an arbitrary number of 
sites. The method enables us to calculate the genuine 
multisite entanglement for moderately large lattices, and 
perform finite-size scaling [17] to predict the measure for 
the infinite two-dimensional square lattice. 

The paper is organized as follows. The succeeding sec- 
tion provides a formal definition of the spin state that 
we consider. Section III contains the theorem proving 
that the spin state is genuinely multiparty entangled. In 
section IV, we present a method to evaluate properties 
of such states on two-dimensional square lattices. In par- 
ticular, a finite-size scaling is presented in Section IV C. 
Wc conclude with some discussions in Section V. 



II. THE SPIN STATE 

An isotropic spin-1/2 system, with dimer coverings, is 
conveniently defined by using a bipartite lattice. A bipar- 
tite lattice is one which is divided into two sublattices, A 
and B, such that each site of sublattice A has only sites 
of sublattice B as its nearest neighbors. Moreover, each 
lattice site is occupied by a spin-1/2 particle. The dimer- 
covering state on such a bipartite lattice can be written 
in the form [18], 

\ip) = ^2 h k( i a, ib)\(ai,h), (a 2 , b 2 )...(a N , b N )) k , (1) 

k 

where each k refers to a single dimer covering ({|a 2 ;, bi)}i) 
on the entire lattice with N spins in each sublattice. 
Here, the "covering function" hk(i a , *&) is a function that 
is isotropic over the lattice sites i a G A,ib 6 B. 

The above states, with short-range dimer coverings are 
often referred as resonating valence bond (RVB) states 
[11, 19]. In quantum spin ladders, these states have been 
argued to be possible ground state approximations of cer- 
tain Heisenberg spin-1/2 systems. Spin systems, such as 
RVB ladders, have been extensively used to study entan- 
glement properties [20, 21]. 
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FIG. 1. (Color online.) A square spin-1/2 bipartite lattice 
with nearest-neighbor dimer coverings. Here, M, M'=6. The 
sublattice A consists of the sites denoted by the lighter cir- 
cles with darker borders, while those in B has darker circles 
with lighter borders. The thick solid lines show the nearest- 
neighbor dimer states (\(o.i,bj))) from a site in sublattice A 
to another in B. The figure represents a possible dimer cov- 
ering. The final state is the superposition of all such dimer 
coverings. 



III. CHARACTERIZING MULTISITE 
ENTANGLEMENT 

A multiparty pure state is said to be genuinely 
multiparty entangled if it is entangled across every pos- 
sible bipartition of the system. Examples of genuinely 
mulitisite entangled states include the Greenberger- 
Horne-Zeilinger [22] and the W [23, 24] states. The 
following theorem shows that the state in Eq. (1) is 
genuinely multisite entangled. 

Theorem: The pure state formed by superpositions of 
dimer coverings, is genuinely multisite entangled, for all 
isotropic spin-1/2 lattices of arbitrary dimensions which 
are periodic or infinite in all directions and all covering 
functions that are isotropic over the lattice. 

Proof. The superposition state of the spin-1/2 lattice 
consisting of 2N particles, in Eq. (1), is a pure state. To 
prove that this superposed state is genuinely multisite 
entangled, we are required to prove that the partial den- 
sity matrix of the state across any bipartition cannot be 
pure. In other words, the density matrix of any p spins, 
formed by tracing the remaining 2N — p spins, is always 
mixed, and hence entangled to the rest. We conveniently 
divide the proof into the two cases where the number of 
spins (a) is finite in at least one part of the bipartition 
and (b) is infinite in both the parts. 

a) Finite case: For a rotationally invariant state, like 
the dimer-covered spin-1/2 state under consideration, it 
is known that the partial density matrix of an arbitrary 
number of spins is also rotationally invariant. Moreover, 
for an odd number of spin-1/2 particles, there is no pure 



quantum state that is rotationally invariant. Hence any 
odd bipartition of the system is always entangled to the 
rest of the system. For example, any single-site density 
matrix is hi, where X is the 2x2 identity matrix, and 
therefore is maximally entangled to the rest of the lattice. 

Let us now consider the case of a bipartition with an 
even (finite) number of spins in one part. Consider any 
set X of an even number of sites. Let the partial density 
matrix of these sites in X, corresponding to the state 
be p( x \ Let us assume that is pure, which would 
imply that is separable, contrary to the statement 
of the theorem. Let X = X' U c, where c contains an 
arbitrary but fixed odd number of sites < \X\. (\S\ de- 
notes the cardinality of the set S.) In particular, \c\ can 
be unity. For an isotropic lattice, we can always find an- 
other equivalent set of spins, Y'Uc, such that \Y'\ = \X'\, 
which is again pure (by the assumption). The strong sub- 
additivity of von Neumann entropy [25] implies 

S(p( x '^) + S(p^) < S(p( x ' uc >) + S(p (y ' uc) ), (2) 

where S(-) denotes the von Neumann entropy of its ar- 
gument. Since pS x ^ and p( Y > are pure, and since von 
Neumann entropy is nonnegative, we have S(p^ x ') + 
S(p( Y )) = 0, which in turn implies that p( x ) and p( Y ) 
are pure. This immediately implies that p^ is pure. This 
is a contradiction as |c| is pure. This part of the proof 
was partially presented in Ref. [20]. 

b) Infinite case: Let us begin with the case of an in- 
finite 2D square lattice partitioned into two half-planes 
by an infinite horizontal line. Let us assume that the re- 
duced states are pure, contrary to the statement of the 
theorem. Let Hp denote the sites in one such half plane. 
Consider now the set Hp U Lh, where Lh is the infinite 
horizontal strip of sites with single-site width, and which 
is directly adjacent to Hp. By isotropy, if p^ HpS > is pure, 
p{HpuL H ) j s a i so p U re. Consequently, we can again use 
strong subadditivity to show that p^ LH ^> is pure. Now if 
p(iff) j s p Urc , 5 hy isotropy, p^v 1 2 ) is also pure, where 
Ly is the infinite vertical strip of sites with single- 
site width, and which is obtained from Lh by rotating 
it by 7r/2 radians around the site with Cartesian coordi- 
nates (pi,j>2)- Again using strong subadditivity, we find 
that p^ 1 ^ 2 ) is pure, which is a contradiction. A similar 
proof works for boundaries which are not straight lines. 

One needs a separate proof for the case when we want 
to prove that the infinite horizontal strip Lh<- having an 
width of r sites is not a product with the remaining por- 
tion of the lattice. Assume, if possible, that p( Ll * T ) is 
pure. Then, by isotropy, p^ L v^ 2 ) is also pure, where 
Lyl' P2 ^ is obtained from Lh*- by a n/2 rotation around 
(pi,P2). Again applying strong subadditivity, we have 
that the reduced state of the 2r spins in L#>- n Lyl^ 
is pure, which is a contradiction (by the (a) part of 
the proof). Similar proofs are possible for other infinite 
strips. 

Just like the proof of the (a) part, the (b) part also 
works, with suitable modifications, for arbitrary isotropic 
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lattices of arbitrary dimensions, e.g. the triangular and 
hexagonal lattices in 2D and the cubic lattice. ■ 



IV. QUANTIFYING MULTISITE 
ENTANGLEMENT 

In the preceding section, we have demonstrated 
that the dimer-covcring superposition on an arbitrary 
isotropic lattice is genuinely multisite entangled. It is in- 
teresting to quantify the amount of genuine multipartite 
entanglement present in such states. This is important 
because the qualitative answer in the preceding section 
indicates only the non- vanishing of the entanglement con- 
tent. However, such an entanglement content can asymp- 
totically graze to zero with increasing system-size. Below 
we show that this is not the case. Specifically, we provide 
an analytical method to calculate genuine multipartite 
entanglement of such dimer-covered superposition states 
and then show that the state under consideration pos- 
sess a relatively high genuine multiparty entanglement 
for arbitrarily large system-size. 

To quantify the genuine multisite entanglement in an 
isotropic spin-1/2 lattice, with dimer coverings, we use 
the generalized geometric measure (GGM) [9]. The GGM 
of an i?-party quantum state is the optimized fidelity dis- 
tance of the state \ipn) from the set of all states that are 
not genuinely multiparty entangled. More specifically, 
the GGM (G(\i/jr))) can be calculated as 

G(\i> R ) = l-A 2 max (\i> R }) (3) 

where A max (|^ fl )) = max \(x\^r)\- \x) is an f?-party 
quantum state with no genuine multisite entanglement. 
The maximization is over all possible |x) states. 

The GGM [9] of an f?-party pure state can be efficiently 
calculated by using the relation 

G(\i> R )) = l-ma,x{\ 2 K . L \KUL = {A u . . . , A R }, KnL = 0}, 

(4) 

where Xk-.l is the maximal Schmidt coefficients in all 
possible bipartite splits K : L of \ip R ,)- 

The computation of the GGM, therefore, depends on 
the efficient generation of arbitrary reduced density ma- 
trices across all possible bipartitions of the spin sys- 
tem. Certain properties like two-point correlations, en- 
ergy, and partition function, of such large superpositions, 
if found to be ground states of a Hamiltonian, can be 
evaluated in some cases by using approximate methods 
like mean field approximations and rcnormalization tech- 
niques [26] . We propose an iterative method, the density 
matrix recursion method (DMRM) for 2D square spin- 
1/2 lattices, with nearest neighbor dimer coverings (i.e. 
hk{ia, ib) = 1 for nearest neighbors, and vanishing other- 
wise) , by which the limitations in computing all possible 
density matrices and its correlation properties of such 
lattice states can be overcome. The method proceeds 
by deriving an algebraic recursion relation between the 



quantum states of small-sized spin systems that can be 
exactly computed. 

We now apply the DMRM to two kinds of 2D systems: 
a) The "perfect" square lattices with an equal even num- 
ber of sites on the horizontal and vertical sides, and b) 
the "imperfect" square lattices having an even number of 
sites on the horizontal side but an odd number on the ver- 
tical one. It is important to note that the total number 
of spin sites is always even. As we see later, the difference 
between the behavior of GGM of perfect and imperfect 
square lattices significantly reduces as the lattice size is 
increased. 



A. Perfect Square Lattice 

We consider a spin-1/2 square lattice with M spins 
along the horizontal and Af' spins in the vertical sides, 
such that M.M' (=2N) is the total number of spins. For 
a perfect lattice with even Af (Af' = Af ), the spin state 
of the system, which we rename as |Af, Af }, can be gener- 
ated by using a recursion of two smaller sized imperfect 
spin states \M-1,M) and \M-2,M). We drop the sec- 
ond Af (=Af' ) from the states, as it remains unchanged. 
For a perfect lattice system, with Af = M' = Af + 2, the 
recursion relation, with open boundary condition, can be 
written as 

\N + 2) = |7V+l)|l} n+2 + |A0|2) n+l,n+2 

= |A0|2) + \M-1)\2) n,n+l |l)n+2, (5) 

where the subscripts correspond to the numbering of the 
sites on the horizontal side on the lattice. For our anal- 
ysis, we consider a spin-1/2 square lattice state, with 
nearest neighbor dimers, that is periodic along the hor- 
izontal axis [27] (cf. [28]). Here, |l)j refers to the state 
corresponding to (Af, M') = (1, Af + 2) at the column i, 
and \2) it j refers to the state with (Af, Af ') = (2, Af + 2) 
at the columns i and j. The state \2)ij is \2)i,j — 
The subtraction in the term [2),*^ removes a repetition in 
the recursion. The recursion can be extended to states 
with the periodic boundary condition as 

\M+2) P = \M+2) hn+2 + \Af) 2,n+l |2)n+2,l- (6) 

The states without subscript P implies non-periodic 
states. The recursion for the corresponding density ma- 
trix is given by 

p ^+ 2 > = p W) + \Af +2)(AA|2,„+i<2U 2 ,i 

+ |AA) 2 ,„ + i|2) 1 ,„ +2 (AA+2||AA)(AA| (2i „ +1) |2)(2| (lin+2) 

(7) 

where p<- A/ '+ 2 *> = \Af + 2) {Af + 2| is the density matrix for 
the non-periodic recursion. Hence using (7), we can gen- 
erate the reduced density matrices across different parti- 
tions that can be used to find the behavior of GGM in 
the system. Tracing over all the spins except the spins 
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at columns n + 1 and n + 2, we obtain the correspond- 
ing reduced density matrix in the case of open boundary 
condition as 

P(n+l.n+2) = ZM 2 )( 2 kn+l,n+2) + ZjV-lPn+1 ® l)(l|(n+2) 
+ (|2)rM-l, ra +2(l|n+2(6vU+l + h.C.), (8) 

where Z^ = (Af\Af), p n +i = tr„[|2)(2| (nj „ +1 )], and 
(6v|n+i = (2\ n , n +i(Af — 1 |A/") . The reduced state in the 
case of periodic boundary condition is 

^(n+l.n+2) = Pn+l.n+2 + tri..„[|^>(^| (2 , n+1) |2)(2| (1 , n+2) 
2,n+l |2)l.n+2 (JV + 2\ + h.c.)} 
~ Vn+l,n+2 "t" ^l(„+l,n+2) T ^2(n+l,n+2) ' ^ ' 

Here. 

^l(n+l,n+2) = ^AT-l|l)(l|(n+l) ® Pn+2^AA-2Pn+l ® Pn+2 
+Zj^- 2 pn+l H j5„+2(|l)(CAA-l|(n+l) ® Pn+2 + fo.C.), (10) 

^2(n+l,n+2) = |2)re+l,n+2(l|»+l(6vU+2 + |2)n+l,»+2 

X y^(J r i|, t+ 2(^-i|»+l + Pn+l ® |l)„+2(CAT|n+l 
i=l 

1 ^ 
+ ^(|J"l)n+l|l)n+2(l| n+ lV(J-i| n +23^-l), (H) 

where yjj = (J\f\J\f — 1)|1). The states 7^) can be recur- 
sively generated, as shown below. To obtain the reduced 
density matrix for a large spin system, the parameters 
that appear upon tracing must be recursively generated 
from the parameters corresponding to smaller systems 
that can be exactly calculated. 

The different "inner products" can be calculated as 
follows. The recursion begins with evaluating (1|2) = 

2i=i a l \ a i}i with {|o!i)}i forming an independent set of 
vectors consisting of certain singlet combinations of an 
(l,A/"+2) spin system, where k has to be found explicitly. 
E.g., |ai) = |l). In general, we can write 

„K-|2)„,„ +1 = (12) 

i 

Using the above relation, we can generate the inner prod- 
uct recursions as follows: 

Z M = (M\N) n 

= Z X Z M ^ + Z' 2 Z M _ 2 + 2(-l)"- 1 

i 

(13) 

>^ = n<ail(l,n-l(A/'-l|JS0 ll „) 



where Z 2 = (2|2) and Aij = (a.i\a.j). Similarly, the other 
states can be generated using subsequent recursions, 

{isr\n+X = (2\n,n+l(N - 1\M) = ^ ■ZjV-i^ln+l, 

i=l 

(15) 

where = „,„+i(2|J*j_i}„ and |7 r ) Il = |1)„. On 

expanding, we can write 

(J~i\n+1 = n.n+1 (2| J r i-l) n = 9k( a k\n+l, (16) 

fc 

where g* = J2 k gjT 1 ^ . Hence all the terms can be re- 
cursively calculated provided the paramters for small spin 
systems can be accurately estimated. The terms a\ = 
1 to A:), of Zi, Z' 2 and Aij need to be exactly calculated. 
The value of A; needs to be determined for a small system 
size by solving the linear equation system (12). 

B. Imperfect Square Lattice 

Let us now consider spin lattices with M sites on the 
horizontal side and M' on the vertical one, where M 
is even, while M' is odd. Here MM' (=2N) and M' 
(=M± 1). For a system containing M = J\f+ 2 spin sites 
along the horizontal side, the recursion of the periodic 
state \Af + 2, J\f + 2 ± 1} p can be written as 

2,n+l|2) n +2,l- (17) 

Here, the constant M' = J\f + 2 ± 1 has been omitted 
from the notation of the state. \2)ij refers to the state 
corresponding to (M, M') = (2, M + 2 ± 1). 

The reduced density matrices can be recursed by using 
Eq. (17). In particular, for obtaining a reduced density 
matrix of the spins at columns (n+ 1, n+ 2), we trace out 
the spins in the columns ranging from 1 to n, and obtain 

Pp\ n+ l, n+ 2) =^l..n[W+2){N+2\ P ] 

= tr 1 .. n [pS/)(A/"| ll „|2)(2| + |A/')(A/'| 2 ,„ + i|2)(2|„ + 2 i i 

l,n|2) n +l,n+2 2,n+l (2|l,n+2 + h.C 

)]■ (18) 
After simplification, the above equation reads 

Pp{n+l,n+2) = ^AA|2)(2|( n+l! „ +2 ) + Zj^-2pn+l ® Pn+2 

+ (|2> n+l,n+2 

(XM\n+l,n+2 + h.C.), (19) 

where the normalization is defined as Zj^ = (Af\Af) = 
Z 2 ^ 2 . The recursion for the state \xj\f) is given by 
(xMn+l,n+2 = (2|i,„ +2 (A/'| 2 ,„ + i|7V)i,„ 
= (2|i,n+2(2|2,3-(2U,„+i|2) 1)2 |2) 3i4 ...|2) n _ 1)n . (20) 

The last recursion can be done numerically to solve for 
(2|i,i+3(2|i+i,i+ 2 |2)i i i + i. 

In the succeeding subsection, finite size calculations 
shows that the behavior of multisite entanglement of the 
two types of lattices merge as the size of the system in- 
creases. 
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FIG. 2. (Color online.) The behavior of GGM in the case of a 
square spin- 1/2 lattice, with increasing total number of spins 
(2N). 



C. Finite size scaling 

We are now ready to calculate the genuine multipar- 
tite entanglement of spin-1/2 square lattices by using the 
DMRM for such lattices. We find the GGM for perfect 
as well as imperfect square lattices and observe that for 
increasing system size, the GGM converges to the value 
0.358. 

In Fig. 2, we plot the GGM with increasing total num- 
ber of spins (2N). Using finite size scaling, the behavior 
of GGM for finite-sized lattices can be used to estimate 
the GGM for an infinite square lattice. The scaling anal- 
ysis gives 



G(|</>}) «G c (|V))±fcN" 



(21) 



where N = 2N and G c (\ip)) is an estimated value of GGM 
for the infinite lattice, based on the average of the last two 
values of GGM given in Fig. 2, with k being a constant. 
The value of estimated by finite-size scaling, using 
G C (|<?W))=0.358, is x =1.69, while fc=0.843. 



V. DISCUSSIONS 

We investigate multipartite entanglement of a large su- 
perposed state consisting of dimer coverings of a spin- 



1/2 isotropic lattice. We first showed that the state, if 
isotropic over the lattice, with either periodic boundary 
conditions or is an infinite lattice, is genuinely multipar- 
tite entangled, regardless of the geometry and dimension 
of the isotropic lattice. To measure its multisite entangle- 
ment content, we have presented a technique, which we 
have referred to as the density matrix recursion method 
for square lattices, to analytically study arbitrary state 
parameters of an arbitrary number of sites, including its 
genuine multisite entanglement. The method was then 
employed to investigate the finite-size scaling behavior 
of the generalized geometric mcrasurc, which is a mea- 
sure of genuine multipartite entanglement, for finite-sized 
square lattices that enables us to to estimate the genuine 
multipartite entanglement of an infinite square lattice. 

Genuine multisite entanglement is potentially a ba- 
sic ingredient in building large-scale quantum comput- 
ers and also in implemntation of multiparty quantum 
communication. The method presented can be a use- 
ful tool if such highly superposed systems are considered 
for performing quantum tasks. Specifically, the iterative 
method can be employed to derive reduced density ma- 
trices that will in turn be fruitful in the calculation of 
nearest-neighbor bipartite entanglement as well as that 
of other two-point correlation functions. 
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